Measurement of the pn^dK + K total cross section close to threshold 
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Measurements of the pd — > p sp dK + K~ reaction, where p Bp is a spectator proton, have been un- 
dertaken at the Cooler Synchrotron COSY- Julich by detecting a fast deuteron in coincidence with 
a K + K~ pair in the ANKE facility. Although the proton beam energy was fixed, the moving target 
neutron allowed values of the non-resonant quasi-free pn — > dK + K~ total cross section to be de- 
duced up to an excess energy e w 100 MeV. Evidence is found for the effects of K~ d and KK final 
state interactions. The comparison of these data with those of pp —* ppK + K~ and pp — » dK + K° 
shows that all the total cross sections are very similar in magnitude. 



PACS numbers: 25.40.Ve, 13.75.Cs, 14.40.Cs 

We have recently published measurements of the dif- 
ferential and total cross sections for the pp — ► ppK + K~ 
reaction at three energies close to threshold A major 
challenge in the analysis was the separation of the con- 
tribution from the production and decay of the meson 
from that of the non-0 component 0]. One of the strik- 
ing features of the non-0 results is the strong attraction 
between the K~ and each of the final protons seen in the 
differential distributions. This also has a major effect on 
the energy dependence of the total cross section, enhanc- 
ing it at low energies. Although tantalizing, these re- 
sults do not, however, resolve the ongoing question as to 
whether the interaction is sufficiently strong to allow the 
K~ to form a bound state with the two protons [H, 0, HI]. 

The isospin dependence of production has been stud- 
ied through an investigation of pd — ► p sp dK + K~ @. By 
identifying the final deuteron and kaon pair and mea- 
suring their momenta, it was possible to construct the 
momentum of the recoil proton p sp to show that it was 
consistent with being a spectator, whose only significant 
participation in a reaction is through a change in the 
kinematics. Interpreting the results in this way, it was 
possible to extract values of the quasi-free pn — ► dK + K~ 
cross section. Moreover, although the experiment was 
carried out at one fixed beam energy, the movement of 
the target neutron enabled data to be obtained over a 
wide range of excess energy e = -Js — rrid — 2ttik on 
an event-by-event basis. Just as in the pp — * ppK + K~ 
case, the shape of the K + K~ invariant mass distribution 



was used to separate the component from the non-0 
background. The invariant K + K~ mass spectrum for 
all events above the threshold is to be found in Ref. [(| 
and from this it is already seen that the non-0 contri- 
bution is a much smaller fraction of the total than in 
the pp — » ppK + K~ case The prime purpose of this 
work is to present the data on the energy dependence of 
the non-0 total cross section up to an excess energy of 
e Ps 100 MeV. 

Unlike production, there are two different 
pn — ► dK + K~ isospin channels. The 1=1 has already 
been investigated in some detail through the measure- 
ment of — > dK + K° at two beam energies, correspond- 
ing to e = 47 and 105 MeV @, U- The identical na- 
ture of the initial protons, combined with angular mo- 
mentum and parity conservation laws, demands that the 
dK + K° final state must contain at least one p-w&ve. At 
low energies this will suppress the 1=1 contribution to 
pn —> dK + K~ compared to / = where there is no such 
constraint. As a consequence, the energy dependence of 
the pn — > dK + K~ total cross section is expected to be 
more complicated than that olpp — > ppK + K~ . 

The investigation was carried using a 2.65 GeV pro- 
ton beam incident on an internal target of the Cooler 
Synchrotron COSY. The experimental details and the 
identification of the dK + K~ candidates were described 
for production Q and so we can here be very brief. 
The forward-going deuteron was measured in the ANKE 
magnetic spectrometer Q and both charged kaons were 
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identified in coincidence on the basis of time-of-flight cri- 
teria. After putting a ±3<r cut around the missing mass of 
the spectator proton, about 4500 p sp dK + K~ events were 
recorded. The background from misidcntified ptt + 7t~ 
events was estimated to be less than 7% and effects from 
this were included in the systematic uncertainties. 

The identification of the residual proton as a specta- 
tor is supported by its momentum distribution shown in 
Fig. [1Ja), which follows well the prediction based upon 
the Bonn wave function llOl . 
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FIG. 1: (a) Spectator momentum distribution of non-0 events 
compared to simulation using the Bonn potential [10]. (b) 
Acceptance-corrected angular dependence of the polar angle 
of the K + K~ system relative to the beam axis in the overall 
cm. system at 12< e <32 MeV. The data are well described 
by 1+0.4 cosfl (solid line). 

The effective target density was determined by measur- 
ing the frequency shift of the stored proton beam as it 
lost energy due to its repeated passages through the tar- 
|. Combining this with a measurement of the 
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beam current, an integrated luminosity of (23 ± 1.4) pb _1 
was found over the 300 hours of data-taking. 

The excess energy determined through the measure- 
ment of the momenta of the deuteron and the two kaons 
had a standard deviation that was typically a c ss 2 MeV, 
which is small compared to the 10 MeV bins that were 
used in the subsequent data analysis. In order to evaluate 
the cross section in one of these e intervals, the geomet- 
rical acceptance, resolution, detector efficiency and kaon 
decay probability were taken into account in a Monte 
Carlo simulation, using the GEANT4 program 13j . The 



fraction of the total luminosity falling within this inter- 
val was estimated from the deuteron Fermi momentum 
distribution predicted using the Bonn potential [Io| . 

In the first step of the analysis, the distributions pre- 
viously published @ were taken as the basis of the sim- 
ulation of the (^-production. For the non-0 component, 
three-body phase space was used. At each excess energy, 
the four independent cm. distributions generated were 
chosen to be the K + K~ invariant mass as well as three 
angular distributions. These were then divided into two 
groups, depending on the value of the K + K~ invariant 
mass, i.e., a 0-rich region where 1.05 < M(K + M~) < 
1.35 GeV/c 2 with the remainder being designated as the 



0-poor region. All distributions were jointly fitted to the 
experimental data and the relative contribution of and 
non-0 production evaluated in order to determine the 
two acceptances. 
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FIG. 2: K K~ invariant mass distributions of the 
pn — * dK + K~ reaction for event samples below and above 
the <f> threshold. The data are shown for excess energy bins 
(a) 42 < e < 52 MeV, and (b) 12 < e < 22 MeV. The non-0 
simulation (dashed curve) includes the effect of the K~d fsi. 
The solid curve includes also a contribution from produc- 
tion, where applicable. 



In the 0-poor region, the polar angle of the K + K~ 
system relative to the beam axis in the overall cm. sys- 
tem shows a forward peak for e < 52 MeV and a typical 
acceptance-corrected distribution is shown in Figfljb). 
The forward/backward asymmetry is evidence for some 
1 = 0/1= 1 interference. The angular distribution for 
the 0-rich region is fairly symmetrical because of the 
1 = dominance and the lack of interference between the 
and the background. Fitting the shape with 1 + a cos 9, 
one finds a = 0.4±0.1 at low energies but a consistent 
with zero at higher e. Taking a to have the polyno- 
mial energy dependence on e, its inclusion increases the 
total acceptance for non-0 production up to 8.5%. In 
addition, as discussed below, the K ± d invariant masses 
deviate from phase space due to the strong final state 
interaction between K~ and deuteron (Fig. [3|). This was 
included through a K~d enhancement factor based on a 
scattering length approximation [8J| . The deviations were 
taken into account iteratively in the simulations in order 
to converge on acceptance-corrected distributions. 

Two typical K + K~ mass spectra from above and be- 
low the threshold are shown in Fig. Rafter making ac- 
ceptance and other corrections. Also illustrated there are 
fits to the and non-0 contributions to the cross section 
where, in the latter distorted three-body phase 

space has been assumed for the dK + K~ final state. In 
general the contribution is well described but the same 
cannot be said for the non-0 distribution. The difficul- 
ties here arise principally from the unknown fraction of 
p-waves that come from the 1 = 1 cross section and the 
influence of the K~d final state interaction. In addition, 
for invariant masses below about 995 MeV/c 2 , the data 
in Fig. HJa) lie well above the simulation. This feature, 
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which is also seen quite clearly in the pp — > ppK + K~ 
data [H, 03], is evidence for a final state interaction in 
the KK subsystem. In the region of the K°K° thresh- 
old, the K + K~ ^ K°K° coupling can lead to a cusp 
effect O. 
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FIG. 3: Ratios of the K ± d invariant mass distributions for 
the same two ranges of excess energy as in Fig. (a) : 42 < 
e < 52 MeV, (b): 12 < e < 22 MeV. The histograms are 
the simulations in the scattering length approximation with 
a = (-1.0 + il.2) fm. 

Of much more importance for the acceptance estima- 
tion is the distortion in the K~d subsystem. Figure [3] 
shows the ratio of the differential cross sections 



RlCd = 



da/dM(K~d) 
da/dM(K+d) 



(1) 



for the non-0 region at the same excess energies as those 
shown in Fig. [2j Here M(K ± d) is the invariant mass of 
the K^d subsystem. Experimental distributions of Rxd 
at both energies show a very strong preference for low val- 
ues of M(Kd), which arises from the interaction between 
K~ and deuteron. A K~d final state interaction factor 
is introduced into the three-body phase space simulation 
of the non-0 contribution by using the scattering length 
approximation 1/(1 — iqa), where q is the K~d relative 
momentum. The (complex) scattering l eng th is believed 
to be of the order of a ps (— 1.0+il.2) fm [16[, which would 
correspond a bound or virtual state with a binding energy 
of e fa 20 MeV. The ANKE pp -> dK+K° data seem to 
be insensitive to the phase of a but they are best fit with 
\a\ ps 1.5 fm [8]. After taking a = (-1.0 + il.2) fm, the in- 
dividual da/dM(K + d) and da/dM(K~d) distributions 
are well described, as is the ratio Rk<i, which is shown 
for the two excess energy intervals by the histograms in 

Fig. m 

In view of the low statistics and the consequent fluc- 
tuations, the non-0 total cross section in a e bin was 
evaluated in two different ways, (a) by subtracting the 
fit to the component in Fig. [2] and summing the re- 
mainder, and (b) by taking the direct fit to the non-0 
part of Fig. O The average of these two values is given 
as the total cross section in Table HI The difference is a 
major contributor to the systematic uncertainties given 
there. 



TABLE I: Total cross section for the non-0 component of the 
pn — > dK + K~ reaction as a function of the excess energy e 
with respect to the dK + K~ threshold. The first error on the 
cross section is statistical and the second systematic whereas 
that on the energy is the bin half-width. The overall ps ±6% 
uncertainty in the luminosity has not been compounded with 
the other errors. 



e 


0"non-</>(tOt) 


(MeV) 


(nb) 


17.1 ±5.0 


1.8 ±0.3 ±0.1 


27.1 ±5.0 


5.9 ± 1.0 ±0.7 


37.1 ±5.0 


12.4 ± 2.5 ± 1.8 


47.1 ±5.0 


16.2 ± 3.6 ± 1.9 


57.1 ±5.0 


27.6 ±4.8 ±3.2 


67.1 ±5.0 


34.9 ±6.7 ±3.4 


77.1 ±5.0 


38.6 ± 10.0 ±5.0 


87.1 ±5.0 


50.2 ± 14.3 ± 7.7 


102.1 ± 10.0 


69.5 ± 17.5 ± 10.7 


The values of the pn - 


-> dK + K~ total cross sec- 



tion presented in Fig. [4] shows a smooth behavior on 
a logarithmic scale. Also shown are results for the 
pp -► ppK+K- P, [13 and pp -> dK+K° 0. 
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FIG. 4: Total cross section for non-0 KK production in 
nucleon-nucleon collisions near threshold. The closed cir- 
cles denote pn — > dK + K~ data (this work) and pp — > dK + K° 
(open triangles Q), whereas the open circles show the re- 
sults for pp^ppK + K~ from Refs. Q, [H The dot- 
ted curve is the best fit of Eq. © to the pp -> dK + K° 
data whereas the solid curve includes also the isospin-zero 
contribution of Eq. (|3]) so as to describe the energy depen- 
dence of the pn — > dK + K~ total cross section. The dashed 
curve represents the subsequent prediction of Eq. (f4]) for the 
pn — > {pn}i = oK + K~ total cross section. 

The isospin dependence of kaon pair production can 
be deduced from 



&(PP —* dK + K°) = (7i, 
a{pn -> dK+ K-) = Ua x 



0o) . 



(2) 
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Interpolating our results to the energies where the 
pp — > dK + K° has been studied 0], we find isospin ratios 
of o-qM = 0.9 ±0.9 at 47 MeV and 0.5 ±0.5 at 105 MeV. 
The large error bars arise from the subtraction implicit 
in Eq. @ and it is hard to draw firm conclusions except 
that (To cannot be much larger than a\ . There might be 
a tendency for a\ to become relatively more important as 
the energy is raised. This is what one would expect from 
the requirement of having a p-wave in the pp — > dK + K° 
final state 0,S|. 

A major influence on the energy dependence of the 
total cross section arises from the K~ d final state inter- 
action which gives the distortion shown in Fig. [3] Al- 
though the closed form description of Ref. [l8( is only 
strictly valid for a real scattering length, to a good ap- 
proximation the energy dependence of the two cross sec- 
tions should be given by 



00 
D 



A e 2 /D , 

A 1 e 3 (D + ie/eo) / D 2 , 

\ 2 

+ y/l + e/e J , 



(3) 



where e ~ 20 MeV. 

The choice of the I — and 1 = 1 coefficients 
A Q = 127 pb/MeV 2 and A x = 1.8 pb/MeV 3 leads to the 
fits to the pp — > dK + K° and pn — > dK + K~ total cross 
sections shown in Fig. 2] The general behavior is repro- 
duced much better than it would be if the K~d final state 
interaction were neglected. 

Another important feature of Fig. fj] is that the 
pp — > ppK + K~ and pn — > dK + K~ total cross sections 
are similar in magnitude. However some allowance has 
to be made for the four-body nature of the ppK + K~ 
phase space. An estimate of this effect can be obtained 
in a simple final state interaction model [191 ]. This pre- 
dicts that 



a(pn 
2 



{pn}i =0 K + K )/ a(pn 
1 - 



dK + K~ 



Ax 
15 



1 

2x 



2y/x~ 



2.i; 



arctan-^/x 



,(4) 



where x = e/B, with B denoting the deuteron binding 
energy. The result of multiplying this ratio by the fit to 
the pn — > dK + K~ total cross section is shown in Fig. 2] 
Using this simple estimate we see that 



cr(pp — > ppK + K )/<r(pn — » {pn}] = oK + K ) w 1.5. 



It is clear from the results presented here that, after 
correcting for the different phase spaces, the total cross 
sections for the pp~>ppK + K~, pn^dK + K~, and 
pp — > dK + K° reactions are very similar in magnitude 
despite the necessity for p-waves in the last case. It 
would be highly desirable to have a common theoretical 
model to describe all three channels. 
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